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Abstract. 

MAGIC is a single-dish Cherenkov telescope located on La Palma (Spain), hence with an optimal 
view on the Northern sky. Sensitive in the 30 GeV - 30 TeV energy band, it is nowadays the only 
ground-based instrument being able to measure high-energy y-rays below 100 GeV. We review the 
most recent experimental results obtained using MAGIC. 
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INTRODUCTION: THE MAGIC TELESCOPE 

The Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope is a last- 
generation instrument for very high energy (VHE) y-ray observation exploiting the 
Imaging Air Cherenkov (lAC) technique. This kind of instrument images the Cherenkov 
light produced in the particle cascade initiated by a y-ray in the atmosphere. Located on 
the Roque de los Muchachos Observatory, in La Palma (Spain), MAGIC incorporates 
a number of technological improvements in its design and achieves the lowest energy 
threshold (55 GeV with the nominal trigger, 25 GeV with the pulsar trigger [1]) among 
instruments of its kind. MAGIC signal digitization utilizes 2 GSample/s Flash Analog- 
to-Digital Converters, and timing parameters are used during the data analysis 112], 
yielding a sensitivity (at a flux peak energy of 280 GeV) of 1 .6% of the Crab Nebula flux 
in 50 hours of observations. The relative energy resolution above 200 GeV is better than 
30%. The angular resolution is ~ 0.1°, while source localization in the sky is provided 
with a precision of ~ 2'. MAGIC is also unique among lAC telescopes by its capability to 
operate under moderate illumination [3] (i.e. moon or twilight). This allows to increase 
the duty cycle by a factor 1.5 and a better sampling of variable sources is possible. 
The construction of a second telescope is now in its final stage and MAGIC will start 
stereoscopic observations in the following weeks. 

MAGIC has been operating since fall 2004, developing a physics program which 
includes both, topics of fundamental physics and astrophysics. In this paper we highlight 
MAGIC latest contributions to ExtragalacticjHJl 0, 0, S @, ES Et [ll] and Galactic 



[H, H H E H [13, H 0, H ImQ astrophysics . 
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EXTRAGALACTIC OBSERVATIONS 



Except for the radio galaxy M 87, all 23 currently known VHE y-ray emitting AGNs 
[|26|l are high-frequency peaked BL Lac object^. Their two-bumped spectral energy 
distributions (SED) are characterized by a second peak at very high y-ray energies. Li 
synchrotron-self-Compton (SSC) models this peak is assumed to be due to the inverse- 
Compton (IC) of electrons, by upscattering previously produced synchrotron photons 
to high energies. In hadronic models, instead, interactions of a highly relativistic jet 
outflow with ambient matter, proton-induced cascades, or synchrotron radiation off 
protons, are the origin of the high-energy photons. Another defining property of blazars 
is the high variability of their emission ranging from radio to y-rays. For VHE y-ray 
blazars, correlations between X-ray and y-ray emission have been found on time scales 
ranging from ~ 10 minutes to days and months (see, e^., Fossati et al. [|27D), although 
the correlations have proven to be rather complicated [281]. 

Here we present selected results from multi-wavelength (MWL) campaigns with 
MAGIC participation and for observations of the AGN Mkn 501, M 87 (February 2008), 
lES 1011-^496, S5 0716-^71, and 3C 279. 



Multi- Wavelength Campaigns 

For an advanced understanding of blazars, coordinated simultaneous MWL observa- 
tions are essential, as they allow the determination of SEDs spanning over 15 orders 
of magnitude in energy. In the 2006/7 season, MAGIC participated in multiwavelength- 
campaigns carried out on the blazars Mkn 421, Mkn 501, PG 1553-1-1 13, lES 1218-1-304, 
IH 1426-1-428, and the radio galaxy M 87. These campaigns involved the X-ray satellites 
INTEGRAL, Suzaku and Swift, the y-ray telescope H.E.S.S., VERITAS, and MAGIC, 
and other optical and radio telescopes. We discuss some of the campaigns in more detail 
here. 

The observations of Mkn 501 in July 2006 revealed the lowest X-ray and VHE state 
ever observed. No variability in VHE y-rays was found, while an overall increase of 
about 50% during one day was seen in X-rays. A one-zone SSC model describes this 
quiescent state of Mkn 501 weU [|29|]. Campaigns on lES 1218-1-304 and IH 1426-1-428 
have been carried out, during both of which significant X-ray variability has been 
observed. Mkn 421 was observed together with INTEGRAL, Suzaku and Swift and 
showed interesting intra-day variability during some of the observation nights. The VHE 
data of these observations are currently being analyzed. A more recent observation of 
a flare of Mkn 421 in summer 2008 could be followed accurately from optical to VHE 
y rays, with the participation of AGILE, GASP-WEBT, VERITAS and MAGIC Q. 
During the first multi- wavelength campaign ever on PG 1553-1-113 in July 2006, its 
VHE emission showed no variability lUj]. lES 1959-1-650 showed VHE data among 
the lowest flux state observed from this object, while at the same time a relatively high 



See |http : / / www .mpp .mpg . de/^rwagner/ sources /[ for an up-to-date source list. 
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FIGURE 1. The night-by-night Ught curve for M87 as measured from 2008 January 30 to 2008 
February 1 1 . Left: The flux in the energy bin 150 — 350 GeV is consistent with a constant emission. Right: 
Flux variations are apparent on variability timescales down to 1 day in the integral flux above 350 GeV, 
instead. The inset shows the light curve above 350 GeV in a 40 min time binning for the night with the 
highest flux (2008 February 1). The vertical arrows represent flux upper limits (95% c.l.) for the nights 
with negative excesses. 



optical and X-ray flux (both Suzaku and Swift) was found [7] . The SED could be modeled 
using similar parameters as needed for the SED measured in 2002, with a slightly more 
compact source and a slightly lower magnetic field. 



Strong Flaring of Messier 87 in February 2008 

The giant radio galaxy M87 has been known as VHE y-ray emitter (Aharonian 
et al. fi321 and refs. therein), and is also one of the best-studied extragalactic black- 
hole systems. To assess variability timescales and the location of the VHE engine 
in M87, the H.E.S.S., VERITAS, and MAGIC coUaborations carried out a shared 



monitoring of M87, resulting in ^ 120 h observations in 2008 1330. Results from the 
entire campaign will appear in a dedicated paper. During the MAGIC observations, 
a strong 8 o signal was found on 2008 February 1 , triggering the other Cherenkov 
telescopes as well as Swift observations. For the first time, MAGIC determined the 
energy spectrum below 250 GeV [|4|], which can be described by a power law with a 
spectral index of F = 2.30 ± 0. 1 Istat ± 0.20syst- We did not measure a high-energy cut- 
off, but found a marginal spectral hardening, which may be interpreted as a similarity to 
other blazars detected at VHE, where such hardening has often been observed [|26|l . Our 
analysis revealed a variable (5.6 a) night-to-night y-ray flux above 350 GeV, while no 
variability was found for 150-350 GeV (Fig. [I]), This fastest variability A? observed so 
far in TeV y-rays in M 87, confirming the E > 730 GeV short-time variability reported 
earlier [|32|l . is on the order of or even below one day, restricting the emission region to 
a size of R< At c5 = 2.6 x 10^^ cm = 2.6 5 Schwarzschild radii (Doppler factor 5), and 
suggests the core of M 87 rather than the brightest known knot in the M 87 jet, HST- 
1, as the origin of the TeV y-ray s. During the MAGIC observations HST-1 was at a 
historically low X-ray flux level, whereas at the same time the core luminosity reached 
a historical maximum (D. Harris, priv. comm.). This strongly supports the core as the 
VHE y-ray emission region. 



The July-2005 Flares of Mkn 501 



MAGIC observed the bright and variable VHE y-ray emitter Mkn 501 during six 
weeks in summer 2005 ^ . In two of the observation nights, the recorded flux (> 4 x that 
of the Crab nebula) revealed rapid changes with doubling times as short as 3 minutes or 
less. For the first time, short 20 min) VHE y-ray flares with a resolved time structure 
could be used for detailed studies of particle acceleration and cooling timescales. 

Interestingly the flares in the two nights behave differently: While the 2005 June 30 
flare is only visible in 250GeV-1.2TeV, the 2005 July 9 flare is apparent in all energy 
bands (120 GeV to >1.2 TeV). Additionally, a photon-by-photon analysis of the July 9 
flare Ig] revealed a time delay between the flare peak at different energies: At a zero- 
delay probability of P = 0.026, a marginal time delay of T/ = (0.030 ± 0.012) sGeV^^ 
towards higher energies was found using two independent analyses, both exploiting 
the full statistical power of the dataset. Several explanations for this delay have been 
considered up to now: (1) Particles inside the emission region moving with constant 
Doppler factor need some time to be accelerated to energies that enable them to produce 
the highest energy 7 rays [5]. (2) The y-ray emission has been captured in the initial 
acceleration phase of the relativistic blob in the jet, which at any point in time radiates up 
to highest 7-ray energies possible [|34|] . (3) A one-zone SSC model, which invokes a brief 
episode of increased particle injection at low energies [35]. (4) An energy-dependent 
speed of photons in vacuum |36|l . as predicted in some models of quantum gravity 
[1371]. When assuming a simultaneous emission of the 7-rays (of different energies) at 
the source, a lower limit of Mqg > 0.21 x lO^^ GeV (95% c.l.) can be established rtgl. 



Blazars Detected during Optical Outbursts 

MAGIC has been performing target of opportunity observations upon high optical 
states of known or potential VHE blazars. Up to now, this strategy has proven very suc- 
cessful with the discovery of Mkn 180 [8], lES 101 1-1-496 [9], and recently S5 0716-1-71 



llllll . The 18.7-h observation of lES 1011-1-496 was triggered by an optical outburst in 
March 2007, resulting in a 6.2 o detection at F>2ooGeV = (1.58±0.32) x lO^^^cm^^s"^ 
[0]. An indication for an optical- VHE correlation is given, in that in spring 2007 the 
VHE 7-ray flux is >40% higher than in spring 2006, where MAGIC observed the 
blazar as part of a systematic search for VHE emission from a sample of X-ray bright 
(^1 kev > 2 /i Jy) HBLs ^ . In April 2008, a high optical state of the blazar S5 07 1 6-1-7 1 , 
triggered MAGIC observations, which resulted in a the detection of a strong 6.8 o sig- 
nal in 2.6 h, corresponding to a flux of F>4ooGeV ~ lO^^^cm^^s"^ flllH . The source was 
also in a high X-ray state LSSJ . The determination of the before-unknown redshifts of 
lES 1011-1-496 {z = 0.21, [ffiand S5 0716-1-71 (z = 0.31, [51) makes these objects the 
third-most and second-most distant TeV blazars after 3C 279, respectively. 
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FIGURE 2. Differential energy spectrum of 3C 279. The grey area includes the combined statistical 
(lc7) and systematic errors. The dotted line shows compatibility of the measured spectrum with a power 
law of photon index a =4.1. The triangles are measurements corrected on the basis of two models for the 
EBL density (see text). 



Detection of the flat-spectrum radio quasar 3C 279 



Observations of 3C 279, the brightest EGRET AGN [|40|], revealed a 5.77 a post- 
trial detection on 2006 February 23 supported by a marginal signal on the preceding 
night The overall probability for a zero-flux lightcurve can be rejected on the 
5.04 a level. Simultaneous optical i?-band observations found 3C 279 in a high optical 
state, a factor of 2 above its long-term baseline flux, but with no indication of short 
time-scale variability. The observed VHE spectrum can be described by a power law 
with a differential photon spectral index of a = 4.1 ±0.7stat ±0.2syst between 75 and 
500 GeV (Fig. [21). The measured integrated flux above 100 GeV on 23 February is 
(5.15±0.82stat±1.5syst) X 10^10 photons cm^^ §-1. 

VHE observations of sources as distant as 3C 279 (z = 0.536) were until recently 
deemed impossible due to the expected strong attenuation of the y-ray flux by the 
extragalactic background light (EBL), resulting in an exponential decrease with energy 
and a cutoff in the y-ray spectrum. The observed VHE spectrum is sensitive to the 
EBL between 0.2 — 2/im. The reconstructed intrinsic spectrum of 3C 279 is difficult 
to reconcile with models predicting high EBL densities (e.g., the fast-evolution model 
of Stecker et al. [41]), while low-level models, e.g. [42!, I43I], are still viable. Assuming a 
maximum intrinsic photon index of a* = 1.5, an upper EBL limit is inferred, leaving a 
small allowed region for the EBL. 

The results support, at higher redshift, the conclusion drawn from earlier measure- 
ments [43] that the observations of the Hubble Space Telescope and Spitzer correctly 
estimate most of the light sources in the Universe. 





FIGURE 3. Left: Skymap of 7-ray candidate events (background- subtracted) around TeV 2032+4132 
for energies above 500 GeV. The MAGIC position is shown with a black cross. Also shown are the last 
positions reported by Whipple and HEGRA. Right: Differential energy spectrum from TeV J2032+4130. 
The shaded area shows the 1 cj error in the fitted energy spectrum. The flux observed by Whipple in 2005 
and in the Milagro scan are marked with squares. The light line shows the HEGRA energy spectrum. 
Theoretical one-zone model predictions are depicted with dashed lines 



GALACTIC OBSERVATIONS 



In this section we review our latest results on Galactic astrophysics. 



The unidentified /-ray source TeV 2032+4130 

The TeV source J2032-I-4130 was the first unidentified VHE y-ray source, and also 
the first discovered extended TeV source, likely to be Galactic [44]. The field of view 
of TeV J2032-I-4130 was observed with MAGIC for 93.7 hours of good-quality data. 



between 2005 and 2007 UM- The source is extended with respect to the MAGIC PSF 
(see Figure |3]). Its intrinsic size assuming a Gaussian profile is CTsrc =5.0±1.7sta±0.6sys 
arcmin. The energy spectrum is well fitted {x^ /n.d.f = 0.3) by the following power 
law: = (4.5 ±0.3) x \0-^^{E/\ re\/)-2 0±0.3 TeV^icm-^s-i. Quoted errors are 

statistical, the systematic error is estimated to be 35% in the flux level and 0.2 in the 
photon index [|l9n . The MAGIC energy spectrum (see Figure [3]) is compatible both in 
flux level and photon index with the one measured by HEGRA, and extends it down to 
400 GeV. We do not find any spectral break, nor any flux variability over the 3 years of 
MAGIC observations. 



Sliell-type Supernova Remnants: Cassiopeia A and IC 433 

We observed the shell-type supernova remnant (SNR) Cassiopeia A during 47 good 
quality hours, and detected a point-like source of VHE y-rays above ~250 GeV [|l4ll 



The measured spectrum is consistent with a power law with a differential flux 

^sys 



at 1 TeV of (1.0±0.1,,flf±0.3™.) x lO^^^ TeV^icm-is^i and a photon index of 




FIGURE 4. Left: Spectrum of Cas A as measured by MAGIC. The upper limits given by Whipple, 
EGRET and CAT are also indicated, as well as the HEGRA detection. The MAGIC and HEGRA spectra 
are shown in the context of the model by |46]. Right: Sky map of 7-ray candidate events (background 
subtracted) in the direction of MAGIC J06 16+225 for an energy threshold of about 150 GeV. Overlayed 
are '^CO emission contours, contours of 20 cm VLA radio data. X-ray contours and 7-ray contours from 
EGRET. The white star denotes the position of the pulsar CXOU J061705. 3+222127. The black dot shows 
the position of the 1720 MHz OH maser. 



r=2.4±0.2 

stat^Q-'^sys- The spectrum measured about 8 years later by MAGIC is consis- 
tent with that measured by HEGRA [|45l] for the energies above 1 TeV, i.e, where they 
overlap (see Figure |4] Left). Our results seem to favor a hadronic scenario for the 7-ray 
production, since a leptonic origin of the TeV emission would require low magnetic 
field intensities, which is in principle difficult to reconcile with the high values required 
to explain the rest of the broad-band spectrum. However, hadronic models [46] predict 
for the 100 GeV - 1 TeV region a harder spectrum than then measured one. 

We have detected a new source of VHE y-rays located close to the Galactic Plane, 
namely MAGIC J06 16-1-225 [15], which is spatially coincident with the SNR IC 443. The 
measured energy spectrum is well fitted (x^/n.d.f = 1.1) by the following power law: 
= (1.0±0.2) X 10-"(£/0.4 reV)-3-i±o.3 TeV-icm-^s-i. MAGIC J06 16-1-225 
is point-like for MAGIC spatial resolution, and appears displaced to the south of the 
center of the SNR shell, and correlated with a molecular cloud [47] and the location of 
maser emission (see Figure|4]Right). There is also an EGRET source centered in the 
shell of the supernova remnant. The observed VHE radiation may be due to ;r°-decays 
from interactions between cosmic rays accelerated in IC 443 and the dense molecular 
cloud. A possible distance of this cloud from IC 443 could explain the steepness of the 
measured VHE y-ray spectrum. 



Wolf-Rayet binaries 

WR stars display some of the strongest sustained winds among galactic objects with 
terminal velocities reaching up to Voo > 1000 — 5000 km/s and also one of the highest 
known mass loss rate M ~ 10^^... 10^^ Mq/ yr. Colliding winds of binary systems 
containing a WR star are considered as potential sites of non-thermal high-energy photon 
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FIGURE 5. Left: Inverse Compton (IC) spectra of WR 147 for orbital phases 0, 0.25, 0.5 and 0.75 fil 
together with MAGIC experimental upper limits. Right: Skymap of y-ray excess events (background 
subtracted) above 150 GeV around Cygnus X-1 corresponding to the flare detected on 2006-09-24. The 
cross shows the best-fit position of the 7-ray source. The position of the X-ray source and radio emitting 
ring-like are marked by the star and contour, respectively. 



production, via leptonic and/or hadronic process after acceleration of primary particles 
in the collision shock (see, e.g., [49]). 

We have selected two objects of this kind, namely WR 147 and WR 146, and observed 
them for 30.3 and 44.5 effective hours, respectively lfl6ll . No evidence for VHE 7-ray 
emission has been detected in either case, and upper limits to the emission of 1.5, 1.4 
and 1.7% (WR 147) and 5.0, 3.5 and 1.2% (WR 146) of the Crab Nebula flux are derived 
for lower energy cuts of 80, 200 and 600 GeV, respectively. These limits are shown in 
Figure [5] (Left) for the case of WR 147, compared with a theoretical model [49,1 . 



Cygnus X-1 is the best established candidate for a stellar mass black-hole (BH) and 
one of the brightest X-ray sources in the sky. We have observed it for 40 hours along 
26 different nights between June and November 2006. Our observations have imposed 
the first limits to the steady 7-ray emission from this object, at the level of 1% of the 
Crab Nebula flux above ~500 GeV. We have also obtained a very strong evidence (4.1 o 
post-trial significance) of a short-lived, intense flaring episode during 24th September 
2006, in coincidence with a historically high flux observed in X-rays [50] and during the 
maximum of the ~ 326d super-orbital modulation [51]. The detected signal is point-like, 
consistent with the position of Cygnus X-1. The nearby radio-nebula produced by the 
jet interaction with the interstellar medium [52] is excluded as the possible origin of an 
eventual putative emission (see Figure [5] Right). 

LS I -1-6 1 303 is a very peculiar binary system containing a main-sequence star together 
with a compact object (neutron star or black hole), which displays periodic emission 
throughout the spectrum from radio to X-ray wavelengths. Observations with MAGIC 
have determined that this object produces 7-ray s up to at least ~4 TeV [|18l] . and that the 
emission is periodically modulated by the orbital motion (/*reV ~ (26.8 ± 0.2) d) [ 24l] 
(see Figure |6] Left). The peak of the emission is found always at orbital phases around 



Compact binaries: Cygnus X-1 and LS I +61 303 
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FIGURE 6. Left: Lomb-Scargle periodogram for LSI +61 303 data (upper panel) and simultaneous 
background data (middle panel). In the lower panel we show the periodograms after subtraction of a 
sinusoidal signal at the orbital period (light line) and a sinusoidal plus a Gaussian wave form (dark line). 
The vertical dashed line corresponds to the orbital frequency. Inset: zoom around the highest peak, which 
corresponds to the orbital frequency (0.0377d^'). Its post-trial probability is ^10^^. The IFS is also 
shown. Right: Pulse profile of the Crab Pulsar above a) 60 GeV and b) 20 GeV, measured by MAGIC. 
The shaded areas show the signal regions for the main pulse (PI) and the inter pulse (P2). 



0.6-0.7. During December 2006 we detected a secondary peak at phase 0.8-0.9. Be- 
tween October-November 2006, we set up a multiwavelenght campaign involving radio 
(VLBA, e-EVN, MERLIN), X-ray (Chandra) and TeV (MAGIC) observations Q. We 
have excluded the existence of large scale (~ 100 mas) persistent radio-jets, found a pos- 
sible hint of temporal correlation between the X-ray and TeV emissions and evidence 
for radio/TeV non-correlation. 



Crab Nebula and Pulsar 

The Crab Nebula is the standard candle for VHE astrophysics and as such, a big 
fraction of MAGIC observation time is devoted to this object. Out of it, we have used 16 
hours of optimal data to measure the energy spectrum between 60 GeV and 8 TeV jl^ . 
The peak of the SED has been measured at an energy E = (77 ± 35) GeV. The VHE 
source is point-like and the position coincides with that of the pulsar. More recently, 
thanks to a special trigger setup, we have detected pulsed emission coming from the Crab 
pulsar above 25 GeV (see Figure [6] Right), with a statistical significane of 6.4 a This 
result has revealed a relatively high energy cutoff, indicating that the emission occurs 
far out in the magnetosphere, hence excluding the polar-cap scenario as a plausible 
explanation for the high-energy origin. This is also the first time that a pulsed y-ray 
emission is detected from a ground-based telescope, and opens the possibility of a 
detailed study of the pulsar's energy cutoff, which will help elucidate the mechanism 
of high energy radiation in these objects. More details can be found at llssll . 

We thank the Instituto de Astrofisica de Canarias for the excellent working conditions 
at the Observatorio del Roque de los Muchachos in La Palma. 
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